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We show that the quahtative picture of the phase diagram which includes the non-uniform 
chiral phase and 2SC superconducting phase is independent of the considered regularization 

'do'- 

Q ' schemes. We also demonstrate that the quantitative results agree with each other reasonably 

o : 

■ for the set of so called "relativistic" regularization schemes. On the other hand the "non- 
O ' 

. relativistic" momentum cut-off is clearly differ from the others. 

Q: 
(N : 

, ■ I. INTRODUCTION 

^: 

fi^ . The phase diagram of strongly interacting matter is in the center of scientific interest for many 
decades. During the last decade there was a great progress in our understanding of this subject which was 



based on the simple observation that Cooper instability and asymptotic freedom jointly lead to the colour 



^ superconductivity phenomenon at asymptotically high densities [1]. However, at moderate densities our 



knowledge is still limited to the model dependent calculations. A variety of phases are possible which 

is a direct consequence of the richness of strong (and indirectly weak) interactions between elementary 
^ ■ 

^ constituents of matter. In particular, the high symmetry groups involved in the description allow the 
00 

O existence of many different types of phases. Of special interest are non-uniform phases. Among those, 

'• 

• ^ ', the LOFF phases of superconductivity [2], Overhauser effect [|3|] or non-uniform chiral condensates [|4|] 
;h are very good examples. 



In this paper, we ask the technical but important question of the regularization dependence of the 
phase diagram of strongly interacting matter which contains the non-uniform chiral condensate. This 
drawback is inferred by the fact that the NIL model describing strong interactions is an effective non- 
renormalizable approximation. Then different regularization schemes lead, in a sense, to different mod- 
els. It is important to check if qualitative results are independent of the regularization scheme, and what 
is their dependence at the quantitative level. 

We consider the phase diagram at finite density which includes chiral uniform and non-uniform 

n 

phases, superconducting 2SC phase and plasma of the free quarks [5]. Herein we checked the phase 
diagram within the Nambu - Jona-Lasinio model against four various schemes, namely 3- and 4-dim 
cut-off, Schwinger and Pauli-Villars regularization. Similar analysis was performed in the case of single 
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non-uniform chiral phase in [|6D. As the main result we confirm that the non-uniform phase exists in all 
considered regularization schemes. 

We expect that the most important differences emerge at zero temperature thus we consider only this 
situation. At higher temperature, the major role would be played by the finite temperature corrections. 
However, inclusion of those is outside the scope of our paper. 



II. MODEL 

The starting point is based on the Nambu - Jona-Lasinio model with two flavours [|70 

H = {^f{if^v+Pjo)^ + G [(iffXK)^ + (WsW)^] + G'(\ff,/Y5X2X'^v) (W5X2?^'^¥c) } (1) 

where V|/ is the quark field, = Cvj/^ is the conjugate field and /j is the quark chemical potential. The 
color, flavor and spinor indices are suppressed. The vector x is the isospin vector of Pauli matrices and 
l/^, A = 2,5,7 are three color antisymmetric SU (3) group generators. The integration = jf dx J (fix, 
where (3 is the inverse temperature and derivative operator 3v = (?3t, V) . Two coupling constants G, G' 
describe interactions which are responsible for the creation of quark-antiquark and quark-quark conden- 
sates respectively. Both couplings are treated as independent. There is also an additional parameter A 
which defines the energy scale below which the effective theory applies. It is introduced through the 
regularization procedure. 

We are working in the mean field approximation within the ansatz [|7|] 



(W) = -^cos^-f, (\i/?Y5X» = -^5«3sin^-f, (\|/X2X'^CY5\|/) = ^5a2 (2) 

which describes three possible phases: the chiral uniform phase {q = 0,M 7^ 0, A = 0), the non-uniform 
chiral phase (^7^0,M7^0,A = 0) and the superconducting phase (^ = 0,M = 0,A7^0). All mentioned 
phases can coexist with each other. Using standard methods one can calculate the thermodynamic po- 
tential iim 

^ M2 |A|2 f fik ^ ^yr f ^ , 
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where the limit of zero temperature was already performed. The last integral is divergent. Before we 
introduce different regularization schemes let us convert equation ^ into another form which much 
better suits our purposes and better underlines the physics of the problem. To reach our goal we translate 
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equation ([3]) into another form 

£^,0 = v/(a'±£o)2 + |A|2, Eo = \JP+M\ (4) 

The first three integrals give finite contributions and only the last two are divergent. Let us note that 
in the absence of superconducting state the next to the last term vanishes and the only divergent con- 
tribution follows from the infinite Dirac sea integral. Additionally the last two integrals are separately 
dependent, the first one on A and the other one on wave vector q. This separation is very convenient for 
the regularization procedure. 



III. REGULARIZATION SCHEMES AND PARAMETERS 

In the first step we expand the last term of equation dH) in powers of the wave vector q. It is well 
known that the parameter at the second order is related to the pion decay constant [|4|] 



/ 7— ^ -12 / -— t£oH (5) 

t±J (271)3 J (271)3 2Ml 

where Mq is constituent quark mass at zero density. The formula for the pion decay constant depends 
on the regularization and is known from the earlier literature (e.g. [18[|). This formula together with 
expressions for the chiral condensate fixes the value of G and A parameters (Table 1). More details 
are discussed in the Appendix A. The coupling constant G' cannot be related to any known physical 
quantity. In the vast literature of color supercoductivity, its value is emplaced somewhere between 
G/2 <G'< 3G/2. 

3D 4D S PV 



A 0.635 1.015 1.086 1.12 
GA2 2.2 3.93 3.78 4.47 
Mo 0.33 0.238 0.2 0.22 



TABLE I: Numerical values of the regulaiization parameter A in GeV and dimensionless quantity GA^ for different 
regularization schemes. In the last row, the value of the quark constituent mass at zero density Mq is given in GeV. 
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Taking into account equation ([5]) one can extract the divergent part of thermodynamic potential dH) in 
the form 

= -M / 7^ {Eto - Eto') - 12 / T^^o. (6) 



We consider four types of different regularization schemes [SD 



{2ny 



3-dim cut-off (3D) A restricts the value of three dimensional momentum. The regularized poten- 
tial takes the form 



(27r)3 V ; J (my 

This regularization was frequently used in previous papers e.g. [|5|]. 



(7) 



• 4-dim cut-off (4D) restricts the value of four-momentum vector in Euclidean space. Using the 
formulae 



(A B)- M'^^ln^O+^' 



we introduce the cut-off parameter to the thermodynamic potential through the equation 

-div 



(271) 



(8) 



where d'^kE = dkod^k. 



Schwinger regularization (S) is based on the formula 

A f°° dz 

In p = / — [exp(x(5 + /£)) -exp(x(A + /£))]. 
B Jo X 

This leads us to the regularized expression for the potential 



div 



d^kE r dx 



(271)4 yi/A2 X 



— < exp 



A \2 



i=± 



.2/||r ?{exp[.(-.^.„^)]}. 



(9) 



This regularization was considered in [|9D where only the single non-uniform chiral phase was 
taken into account. 

Pauli-Villars regularization (PV) introduces an arbitrary number of coupling constants Q and 
mass regulators Ma combine in such a way that divergent potential ^ becomes finite. In the 
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first step we regularized potential by the 3 -dim cut-off and then expand the result around the large 
value of the parameter A. The values of coupling constants Ca , Ma are already set by the conditions 
which follows from the calculation of the chiral condensate (vfA[r) and pion decay constant (the 
appendix). Then final expression for the thermodynamic potential in this scheme reads 

3 ^ 4 m2 1 /|A|4 , ,7 7\ ^ m2-|-|A|2 



, + |A|2 

|A|2 3 M2 + |A|2 1 3 4/ M2+|A|2 

+ yiEC-^'S'^JTrT^ + S^iiEcX^ln^^-ln-^), (10) 

where Ca^Ma are given by equations (fT3l) in the appendix. 



IV. RESULTS 



We minimize the thermodynamic potential 



4G 4G' 2M^ 




T7A=0 



(11) 

with respect to mass M, wave vector \q\ and gap parameter A as a function of chemical potential. The 
last term of (fTTI) depends on the scheme as was described in the previous section. Appropriate formule 
for Q.'^ll'"'^ are given by equations (HElllIIOl). 

As already discussed G,A coupling constants are determined by the values of the pion decay con- 
stant and the size of the chiral condensate at zero density (Table I). The remaining constant G' is es- 
sentially unknown. For the presentation of the result we assume G' = 0.75G which follows from the 



Fierz transformation lIlOll . The choice of another ratio G'/G does not influence our analysis of the result 
dependence on regularization schemes. The values of constituent mass, wave vector and gap parameter 
as a function of chemical potential for different regularization schemes are given in figures 1 - 3. As can 
be seen in all regularization schemes, there is the same pattern of the phase transitions. From uniform 
to non-uniform chiral phase and then to superconducting phase, all transitions are first order and exis- 
tence of non-uniform phase is then independent of the considered regularization schemes. However, the 
strengths of the transition depend on the regularization scheme. This is particularly visible in Fig. 3. 

Some quantitative features change with the chosen scheme. One can find that schemes cluster in 
two groups which one can call "relativistic" schemes (4D, S, PV) and 3D cut-off. However, let us note 
that the distinction between relativistic and non-relativistic schemes has no deep meaning because the 
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FIG. 1 : The constituent mass M as a function of chemical potential fi in different regularization schemes in MeV 
units. 

thermodynamic system singles out one reference frame. Comparing the quantitative results, we consider 
values of constituent mass, wave vector, gap parameter, the strengths of the first order phase transitions 
and critical chemical potential. The position of the transition from uniform chiral to nonuniform chiral 
phase is the most resistant against the choice of the regularization scheme. First critical potential changes 
within the range of 5 per cent for relativistic schemes, while in the case of 3D within 18 per cent. 
The position of the second transition, for relativistic schemes changes within the range of 19 per cent, 
including 3D cut-off within 34 per cent. For 4D, S, PV schemes the range of variability of the parameter 
M value is about 20 per cent at /u equal to zero, and about 60 per cent at chemical potentials corresponding 
with first transition. The range of variability of parameter q value is 50 per cent at chemical potentials 
corresponding with first transition. Comparing M and q at second transition is unsuitable because of 
changeability of second critical potential. The dependence of the parameters M and q on the chemical 
potential is the same for the different regularizations. With increasing chemical potential, the value of q 
is growing, the M value is declining. Relatively the least correlated is the dependence of the parameter 
A on the regularization scheme. However, again the value of the gap A is increasing with /j, independent 
of the regularization choice. The values of the critical chemical potentials of the phase transitions are 
given in Table II. 
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FIG. 2: The value of the wave vector ^ as a function of the chemical potential /u in the different regularization 
schemes in MeV units. The values of /n for non-zero q describe the range of the existence of the non-uniform 
chiral phase. 

fi, 3D 4D S PV 
Ch/NCH 0.311 0.274 0.263 0.268 
NCH/2SC 0.373 0.330 0.296 0.281 

TABLE II: Numerical values of the critical chemical potentials jUc in MeV for uniform and non-uniform chiral 
phase transition (Ch/NCH) and non-uniform chiral and superconducting phase transition (NCH/2SC) for different 
regularization schemes. In 3D, 4D, S schemes the 2SC phase appears already at Ch/NCh phase transition and both 
phases coexist. However, the value of the gap parameter is relatively low in the coexistence region. In the PV 
scheme, 2SC phase appears only at the NCh/2SC phase transition. 



Strengths of phase transitions depend on the regularization scheme. The strongest phase transitions 
are in 3D cut-off, the weakest in the Schwinger proper time regularization. In the case of the transition 
to the 2SC phase the jump of the gap ranges from 28 MeV for 3D cut-off to 1 MeV for Schwinger 
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regularization. There is still a possibility of the coexistence between chiral and superconducting phase. 
It occurs in all schemes with the exception of Pauli - Villars. Any conclusion which follows from this 
phenomenon is thus model dependent. Finally, let us note that the phase diagram depends also on the 
value of G' . Its influence is the same for all regularization schemes. The larger value of G' , the shorter 
the range of non-uniform chiral phase. 

This behavior is understandable because larger G' strengthens diquark interaction which dominates 
over quark - antiquark interaction. Only in the Pauli - Villars scheme, the non-uniform chiral phase 
vanishes for G' = 0.83G, and the phase transition to superconducting phase at ^ = 0.26 GeV takes place 
directly from the uniform chiral phase.However, this value of the critical chemical potential is rather low 
which questions the physical sensibility to set G' = 0.83G in PV scheme. 



V. CONCLUSIONS 



We perform the analysis of the phase diagram of strongly interacting matter in the Nambu - Jona- 
Lasinio model which includes non-uniform chiral phase and superconducting 2SC phase in different 
regularization schemes. We confirm that the qualitative features of the phase diagram is independent 
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FIG. 3: The gap parameter A as a function of chemical potential fi in different regularization schemes in MeV 
units. 
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FIG. 4: The phase diagram at one plot in 4D cut-off regularization schemes (in MeV units). 

of considered regularization schemes. The generic phase diagram is shown in Fig. 4 (in 4D cut-off 
regularization). 

The quantitative results (values of constituent mass, wave vector, critical chemical potential) match 
satisfactorily within "relativistic" schemes. Only results obtained with 3D cut-off differ widely from 
the others. Unfortunately, there is no general argument which scheme better suits the task of the phase 
diagram analysis. Neither the relativistic approach (thermodynamic systems single out prefered refer- 
ence frame) nor gauge independence (the NJL model is not a gauge theory) favor any scheme in the 
present considerations. Additionally, the gap parameter A as well as the magnitude of the jump in the 
gap parameter shows the clear dependence on the choice of the regularization scheme. The differences 
appear not only between the "relativistic" and 3D cut-off schemes but also within the set of "relativistic" 
regularizations. These findings tell us that one can set the magnitude of the gap parameter in the MeV 
scale but its precise value in the large extand is an unknown quantity. From the other hand the general 
qualitative pattern that the gap parameter increases with increasing value of the chemical potential /j is 
independent on the choice of the regularization scheme. 

The size of the non-uniform phase depends on the relative strength of G' and G coupling constants. 
The larger G' constant the shorter range of the non-uniform phase. This conclusion is also independent 
of the regularization scheme. 

Finally, we find that in the Pauli - Villars scheme, in contrast to the other schemes, there is no coexis- 
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tence region of the non-uniform and 2SC phases. Thus such a coexistence remains an open question. Let 
us stress at the end that our analysis does not prove that the non-uniform chiral phase exists. However, it 
shows that the main features of the phase diagram which includes non-uniform phase are robust against 
the choice of the regularization schemes. 



VI. APPENDIX A 



Two parameters of the NJL model (G, A) are fixed by two physical quantities: the pion decay constant 
f-j^ = 93 MeV, and the quark condensate density (uu) = (dd) =-(250 MeV)^. These quantities are 
functions of Mq and A, and can be calculated in the framework of the NJL model, as was done in [|8|]. 
Alternatively one can use the decay constant gp for the p — > 27i process instead of the quark condensate 
value, proposed in [[llj]. Now using the self-consistency condition, Mq = —4G{uu), that links G and A 
with Mq, we get values of G and A. 

In different regularization schemes, one has 
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a=l 
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21 



M2 



a=l 



M^ 
/"a 



with three conditions impose on regularization parameters 

3 3 3 

l+£Ca = 0, M2+£CaM2=0, M"" + CaM^^ = 0. 
a=l a=l a=l 



(12) 



Conditions (1121) are solved by the formulae: 



Ci = -3, C2 = 3, C3 = -l, =m2+a2, =m2 + 2A2, M^=M^ + 3A^ (13) 
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The values of G, A for each regularization scheme are given in Table I. 
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